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ABSTRACT 

This paper continues a previous study of neutron stars with positive polar-cap 
corotational charge density in which free emission of ions maintains the surface electric- 
field boundary condition E • B = 0. The composition of the accelerated plasma on any 
subset of open magnetic flux-lines above the polar cap alternates between two states; 
either protons or positrons and ions, of which the proton state cannot support electron- 
positron pair creation at higher altitudes. The two states coexist at any instant of time 
above different moving elements of area on the polar cap and provide a physically 
consistent basis for a description of pulse nulls and subpulse drift. In the latter case, it 
is shown that the band separation P3 is determined not by the E x B drift velocity, as is 
generally assumed, but by the diffusion time for protons produced in reverse-electron 
showers to reach the region of the atmosphere from which they are accelerated. An 
initial comparison is made with the survey of subpulse drift published by Weltevrede 
et al. 
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1 INTRODUCTION 

It is now generally assumed that electron-positron pair cre- 
ation immediately above the magnetic poles is essential 
for coherent radio emission in pulsars. Muslimov & Tsy- 
gan (1992) have made an important contribution in their 
recognition of the significance of the Lense-Thirring effect 
in the generation of adequately strong electric fields (see 
also the paper of Muslimov & Harding 1997). It is remark- 
able, and must be almost unique, that a general-relativistic 
effect changes the order of magnitude of a field component 
in a system that can be described, ostensibly, perfectly well 
in Euclidean space as in the paper of Goldreich & Julian 
(1969). 

At least initially, finding the acceleration field Ey which 
is the component locally parallel with the magnetic fiux 
density B, can be regarded as a problem in electrostatics 
within a volume in the corotating frame of reference hav- 
ing well-defined boundary conditions. This is defined by the 
neutron-star surface and the surface separating open from 
closed magnetic flux lines and is assumed to have the elec- 
tric potential boundary condition "I> = 0. It is possible that 
this surface is actually a current sheet of some complexity 
(see, for example, Arons 2010), but we shall assume that, 
near the neutron-star surface, the time-dependence of Ey 



* E-mail: p.jonesl@physics.ox.ac.uk 



within it is not sufficiently rapid to hinder the charge move- 
ment necessary to maintain the boundary condition. Thus 
there are three possible cases depending, firstly, on the re- 
lation between rotation spin 17 and B and hence the sign of 
the polar-cap Goldreich- Julian charge density po, and sec- 
ondly on the neutron-star surface value of Ey. They are: 

(i) n ■ B > 0, po < 0, E|| = with electron acceleration; 

(ii) n • B < 0, po > 0, E|| =0 with positron and positive 
baryon acceleration or (iii) f2 ■ B < 0, po > and Ey 7^ 0. 
Throughout this paper, for brevity, these conditions will be 
referred to as cases (i), (ii) or (iii). There is no further pos- 
sible case because free emission of electrons is always possi- 
ble at neutron-star surface temperatures. A previous paper 
(Jones 2010a, hereafter Paper I) reported some preliminary 
work on the relation between these sets of boundary condi- 
tions and observed phenomena, in particular, the existence 
of pulse nulls in a significant fraction of radio pulsars. There 
seems to be no reason why neutron stars satisfying (i), (ii) 
and even exceptionally (iii) should not exist and it is the 
purpose of the present paper to continue the attempt to see 
if features of the predicted plasma acceleration in these cases 
correlate with observed properties of subsets of the pulsar 
population. 

There can be little doubt that in this context nulls and 
subpulse drift are important phenomena. Although there 
must be reservations about assigning undue weight to the 
properties of a single neutron star, nulls in the isolated pul- 
sar PSR B1931-1-24 are informative. Kramer et al (2006) 
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found that spin-down in the on-state is approximately twice 
as fast as in the off-state. It is difficult to see, in a pulsar 
of this age (> 10^ yr), how the geometrical shape of the 
acceleration volume or the magnitude of the current density 
J within it can change by a factor of this order in less than 
the rotation period P. Time-varying fields near the light 
cylinder can further accelerate ultra-relativistic particles of 
both signs and therefore the most obvious explanation for 
the change in spin-down torque is that the particle and field 
components of the magnctosphcric momentum density and 
stress tensor near the light cylinder differ between on and 
off-states. A cessation of pair creation during nulls would 
greatly change the charged particle number density even 
though J remains essentially unchanged. In this connection, 
it is interesting that a recent paper by Lync ct al (2010) 
suggests quasi-periodic switching between two different spin- 
down rates as the origin of the peculiar timing residuals seen 
in many pulsars. 

The presence of nulls might be thought of as evidence 
that a pulsar is in the final phase of evolution prior to com- 
plete cessation of radio emission. However, we suggest that 
this view is not consistent with a specific feature of the 63 
rmlling pulsars listed in Tables 1 and 2 of the paper by Wang, 
Manchester & Johnston (2007). The maximum potential dif- 
ference ^max available for acceleration at the polar cap is 
proportional to BdP~^, in which Bd is the effective dipole 
field inferred from the spin-down rate. It is approximately 
independent of the ratio of Bd to the actual polar-cap field 
B. The distribution of this quantity for all radio pulsars in 
the ATNF catalogue (Manchester et al 2005) has a relatively 
sharp cut-off at 2.2 x 10^^ G s^^, equivalent to the existence 
of a well-defined death line in the distribution of Bd versus 
P. It might be expected that nulls should be seen only as a 
pulsar rotation slows so that it approaches this value. But 
the form of the distribution of BdP~^ for the nulling pulsars 
listed by Wang et al is broadly indistinguishable from that 
of the whole ATNF catalogue and instead is consistent with 
nulls being a long-term property of a certain sub-set of radio 
pulsars. It is, of course just possible that the observed sharp 
cut-off is itself a statistical fluctuation and that individual 
pulsar deaths actually occur at values of BdP^^ through- 
out the whole distribution because some other variable is 
involved, such as fiux-line curvature. But we shall discount 
this possibility and in the present paper adopt the view that 
pulsars with boundary condition (ii) and, exceptionally, (iii) 
form the sub-set that null. 

A recent large-scale survey by Weltevrede, Edwards & 
Stappers (2006) has revealed that subpulse drift is a common 
phenomenon. Under the assumption that electron-positron 
pairs are the source of the coherent radio emission, it im- 
plies that compact zones of pair creation exist and move in 
an organized way within the open magnetosphere area of the 
polar cap. This motion has been assumed to be an E x B 
drift velocity under the case (iii) surface electric field bound- 
ary condition, following the original paper of Rudcrman & 
Sutherland (1975). There have been many later papers on 
this subject and we refer to Gil, Melikidze & Geppert (2003) 
for recent developments which have sought to refine the cal- 
culation of the drift velocity. The problem is that the case 
(iii) boundary condition is implausible as a common prop- 
erty of the pulsar population and that cases (i) and (ii) have 
not hitherto provided any immediate and obvious basis for 



the formation of organized subpulses. It was shown in Pa- 
per I that short time-scale instability in the composition of 
accelerated plasma exists in case (ii). In this paper, it is 
shown to be a plausible mechanism for subpulse formation 
and drift. 

Under the assumption of an actual dipole field, pair cre- 
ation is possible, in principle, in all observed pulsars through 
the inverse Compton scattering (ICS) of polar-cap photons 
by accelerated electrons or positrons (Hibschman & Arons 
2001, also see Fig. 1 of Harding & Muslimov 2002), but cur- 
vature radiation (CR) can produce pairs only in those with 
high values of BdP~^. However, there is a problem, noted 
by Harding & Muslimov, in the formation of a dense pair 
plasma because the number of ICS pairs formed per primary 
electron or positron accelerated can be smaller than unity 
(see also Fig. 8 of Hibschman & Arons). They suggest that 
the pair density required for coherent radio emission may be 
far smaller than previously thought. Although higher mul- 
tipole field components may be present in most pulsars, so 
increasing flux-line curvature, the existence of this problem 
must remain a matter of concern perhaps even leading to 
doubts about the role of pair plasma as the source of coher- 
ent radio emission. 

The existence of solutions of the electrostatic problem 
under boundary conditions (ii) or (iii) is no more than a 
preliminary because we are concerned in this paper with the 
presence of instabilities which arise from the reverse electron 
flow at the polar cap. It might be questioned whether or 
not the charge density on the surface separating open from 
closed magnetic flux-lines needed for the condition $ = 
can be maintained in the presence of instability. But we 
shall find that, even at short time-scales, instability princi- 
pally changes the nature of the particles accelerated rather 
than the current density J and the acceleration field. The 
instabilities considered here are not, of course, a feature of 
case (i) in which electrons are the primary component of J. 
Establishing instability in this case is a purely electromag- 
netic problem which has been considered by Levinson ct al 
(2005), Melrose & Luo (2009) and Reville & Kirk (2010). 
Non-stationary flow in case (iii) has been investigated by 
Timokhin (2010). 

The properties of the condensed matter surface are nat- 
urally important in cases (ii) and (iii), in particular, the 
production of protons by the reverse flux of electrons de- 
scribed previously (Paper I; see also Jones 1981). This is a 
characteristic of electromagnetic showers that is usually of 
little practical significance. The form of the shower depends 
principally on electron-photon interactions but shower pho- 
tons also interact directly with nuclei. The cross-section is a 
maximum with the excitation of the giant dipole resonance 
(GDR), a broad collective state, at a photon momentum 
A: w 40 mc. (Electroproduction cross-sections are smaller by 
a factor of the order of the fine-structure constant and can 
be neglected.) Photon track length per unit interval of k in 
a high-energy shower is oc k^^ so that to a good approxima- 
tion, photoproduction of baryons can be assumed to occur 
entirely through GDR formation. Known cross-sections and 
electromagnetic shower theory allow the calculation of pro- 
duction rates, and we define Wp as the number of protons 
produced per unit primary shower energy. Protons are ini- 
tially of the order of a few MeV but are very rapidly moder- 
ated to thermal energies without further strong interaction. 
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and then diffuse to the surface with a time delay that is of 
prime significance for the stabihty of plasma acceleration. 
For further details we refer to Paper I, also to Jones (2010b) 
for the cross-sections at high values of B for processes of sec- 
ond order in electron-photon coupling. In the early stages of 
acceleration, partially-ionized atoms interact with the polar- 
cap blackbody radiation field and this is the source of the 
reverse-electron flux that is considered here. It has proved 
convenient to combine rates for this process with values of 
Wp by defining, for a particular pulsar, the parameter K 
which is equal to the number of protons produced per unit 
nuclear charge accelerated. 

Instabilities in solutions of the timc-indcpcndcnt elec- 
trostatic problem referred to above exist as a consequence 
of proton formation and we propose that they are the basis 
for the nulling phenomenon and for subpulse formation and 
drift. The complexity that arises is unfortunate because it 
limits what can be achieved in terms of a physical theory 
of the acceleration process. Thus we shall be able to ex- 
pose the general properties attaching to cases (ii) and (iii) 
but are not always able to give quantitative predictions. We 
shall show in Section 2 that general consideration of polar- 
cap parameters rule out the possibility that the relatively 
numerous subsets of pulsars that exhibit nulls and subpulse 
drift belong to case (iii). Therefore, sections 3 and 4 of this 
paper are restricted to further consideration of case (ii) for 
which short time-scale instability in plasma acceleration has 
been described previously in Paper I. Its properties arc sum- 
marized in Section 3. The previous treatment of medium 
time-scales was, however, inadequate and contained an er- 
ror. The appropriate analysis is given here in Section 4. The 
relations between these instabilities and the observed prop- 
erties of nulls and of subpulse drift arc described in Section 
5 and, in particular, the relation beween the subpulse band 
separation periodicity Ps and proton diffusion time is given. 



2 POLAR-CAP PARAMETERS 

Many of the polar-cap properties and parameters that will 
be required for the arguments of Sections 3-5 have been given 
previously in Paper I. The properties of the polar cap atmo- 
sphere are of particular importance and will be further con- 
sidered in Section 3. But it will be convenient to summarize 
the remainder here with some additions. 

The most basic parameter is the ion number density 
A'^ of the condensed matter at zero pressure. The magnetic 
dipolc fields Bd inferred from pulsar spin-down rates vary 
over about five orders of magnitude, but the median value 
for the 63 nulling pulsars listed by Wang et al is 2.8 x 10^^ 
G, which is significantly larger than for the whole ATNF 
catalogue. It is also possible that the actual polar-cap field B 
is larger than Bd. For these reasons, we adopt the expression. 



TV = 2.6 X 10^^'Z""-^b!-^ 



(1) 



fitted to values computed by Medin & Lai (2006), primarily 
for Bi2 > 10, where B12 is the magnetic flux density in units 
of 10^^ G and Z is the atomic number. The convenient unit 
of depth below the surface at z = is the radiation length. 
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defined here in terms of the zero-field Bcthc-Heitler 
bremsstrahlung cross-section with screening factor modified 
for the neutron-star surface density (see Paper I). 

The critical temperature above which the ion ther- 
mal emission rate is high enough to maintain the case (ii) 
boundary condition is related to the cohesive energy Ec by 
ksTc w 0.025£c (see the discussion of atmospheric prop- 
erties in Section 3). Cohesive energies have been calcu- 
lated by Medin & Lai as functions of B. For Z = 26 and 
B12 ~ 10, their value is in good agreement with that ob- 
tained by Jones (1985) using a different representation of 
the three-dimensional condensed matter state. In the inter- 
val 10 < B12 < 100, their values can be fitted by the expres- 
sions Ec = 0.016Bh^ keV for ^ = 6 and Ec = O.WBli^ keV 
for Z = 26 giving, 

Tc = 4.6 X 10^5^2^ K, Z = 6 

= 4.6 X W*B%^ K, Z = 26. (3) 

These are to be seen in relation to other temperatures that 
are significant. Paper I contained a calculation of the rate of 

proton formation in the electromagnetic showers formed by 
reverse electrons incident on the polar cap. From the energy 
flux needed to produce a proton current density = poc, we 
can infer a maximum steady-state polar-cap temperature, 

1/4 



(-BcosV)(l - t) 



K. 
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In this expression, Wp is the number of protons produced per 
unit primary shower energy. We can approximate it initially 
by Wjf ^ which was obtained in Paper I by using the zero- 
field Bethe-Heitler pair creation cross-section with screening 
modified for the condensed matter density at the neutron- 
star surface. Its values, given in Table 1 of that paper, can 
be summarized conveniently in the intervals B12 = 10^ — 10^ 
and Z = 10 — 26 by the expression 

(5) 



= 3.9 X W-*{Z);rB'l2 {mc 



in which the nuclear charge is the average in the region 

of the shower maximum. The angle tp is that between f2 
and B. The general- relativistic correction contained in the 
surface value of the Goldreich- Julian charge density po is 
K (see Muslimov & Harding 1997), and a is Stefan's con- 
stant. Equation (4) also contains a further temperature, 
Tres, which is the polar-cap temperature in the absence of 
any reverse electron or photon energy flux (approximately 
the observed whole-surface temperature corrected to the lo- 
cal proper frame). The presence of Tres assumes that there 
is a constant heat flow to the surface driven by the very 
much higher temperature of the inner crust. With Tres = 
and K = 0.15, we find, 



T^„, =5.1 X 10'(Z)°^'b! 



,a2 ^ -cosV> y/^ 



K 



(6) 



Equating this with Tc allows us to estimate the minimum 
polar-cap magnetic fiux densities necessary to sustain the 
case (iii) boundary condition E|| ^ 0. These are, 



181 



327 



-cosVX'''* 



P 

■ cos V"^ 



(2) 



(^) Z = 26. (7) 

Comparison with the median value of Bdi2 = 2.8 for the 
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Table 1. The table gives values of Wp for high magnetic flux 
densities B ^ Be and nuclear charge Z. The effective total cross- 
sections for pair creation in the GDR photon energy band 
are for transverse momenta below the lowest magnetic conversion 

threshold at k± = 2 mc and arc in units of barns. They are 
estimates obtained from Table 3 of Jones (2010b) by averaging 
over photon polarization and over photon transverse momenta 
k± = 1.0 and 1.5 mc, and are equivalent to a mean free path for 
Coulomb pair creation lengthened by a factor rjp compared with 
that for the modified zero-field Bethe-Heitler cross-section. The 
Landau-Pomeranchuk-Migdal effect is not significant in the GDR 
photon energy }>aiid 



z 


BB-^ 


bn 


rip 


(mc^)"-'- 


10 


1 


0.27 


1.14 


2.3 X 10"* 




3 


0.071 


2.95 


3.4 X 10^'' 




10 


0.017 


6.3 


3.6 X 10^"' 


18 


1 


0.89 


1.35 


1.8 X 10""' 




3 


0.23 


3.8 


3.2 X 10"'' 




10 


0.055 


9.7 


3.8 X 10-4 


26 


1 


1.86 


1.45 


1.3 X 10"* 




3 


0.48 


4.2 


2.7 X 10--* 




10 


0.11 


11.3 


3.4 X 10"" 



63 nulling pulsars listed by Wang ct al shows that case (iii) 
can be widely realized only with implausibly large deviations 
from a central dipole field. 

However, the screening-modified zero-field Bethe- 
Heitler pair creation cross-section is not obviously valid 
at magnetic flux densities of the order of the critical field 
Be = 4.41 X 10^^ G. Thus we have been obliged to cal- 
culate the second-order bremsstrahlung and pair creation 
cross-sections using Landau function solutions of the Dirac 
equation. The photoproduction of protons by giant-dipole 
resonance (GDR) decay is determined by the total photon 
track length in the GDR band, centred on a momentum 
fc = 40 mc, which occurs almost entirely in the late stages 
of shower development. The track length at these energies 
is limited by Coulomb and magnetic pair creation, also by 
Compton scattering the effect of which is almost always to 
scatter the photon so that its transverse momentum compo- 
nent k± (perpendicular to the field) exceeds the threshold 
for magnetic conversion to electron-positron pairs. We refer 
to Sections 2 and 3 of Paper I for a more detailed description 
of these processes. Approximate values of the Coulomb pair 
creation cross-section at magnetic fields B ^ Be are given in 
Table 3 of Jones (2010b) and are the basis for the values of 
Wp given here in Tabic 1. These arc not easily summarized as 
a simple expression analogous with equation (5). The cross- 
section a,t B = lOBc is at least an order of magnitude smaller 
the the modified zero-field Bethe-Heitler cross-section but 
the effect on Wp is not large because the photon track length 
in this region is limited by Compton scattering. Substitution 
into equation (4) gives values of Tmax that typically are re- 
duced by a factor of approximately 0.9 from those of equsr 
tion (6). But the complexity of the second-order processes 
at B ^ Be is such that we have not reconsidered shower 
development and, specifically, have not allowed for the pro- 
duction by cyclotron emission or Coulomb bremsstrahlung 
of photons with k± above the thresholds for magnetic con- 
version. This becomes significant at B = lOBe, as shown 



in Tables 2 and 4 of Jones (2010b), and the high magnetic 
conversion transition rates reduce GDR-band photon track 
length in the shower. The reduced Wp values increase Tmax- 
Owing to this complexity at B ^ Be there are inevitable 
uncertainties in our estimates of Wp, but we believe that 
they do not seriously invalidate the estimates of the min- 
imum polar-cap magnetic flux densities needed to support 
the case (iii) boundary condition and our conclusion that it 
can exist only in a very small subset of pulsars. Our conclu- 
sions drawn from equation (3) arc also independent of the 
spectrum of the reverse electron-photon flux because Wp is 
almost completely independent of primary shower energy 
provided that is large compared with the GDR energy. 

There appear to be no published estimates of the melt- 
ing temperature of condensed matter that are speciflc to 
very high magnetic fields. Consequently, we are obliged to 
adopt the standard one-component Coulomb plasma expres- 
sion (see, for example, Slattery, Doolen & DeWitt 1980) 
which, with equation (1), gives, 

Tm = 1.0 X 10*Z^Z-°-2^B?2* K, (8) 

in terms of an efi^ective valence charge Z„ . This latter param- 
eter represents the fact that the deeply-bound Landau states 
certainly do not participate significantly in the melting tran- 
sition, but its estimation at higher values of Z is quite dif- 
ficult. It is possible that the work of Potekhin, Chabrier & 
Yakovlev (1997; see Fig. 1) could provide guidance although 
it is at zero field and was directed toward a different prob- 
lem. On this basis, we assume for Z = 26 a value in the 
interval = 10 — 15. In a typical polar-cap field, B12 = 10, 
the melting temperature is as low as Tm = 6 x 10^ K for 
Z„ = Z = 6 but exceeds 10** K for Z = 26. Thus the state 
of the polar cap may be a sequence of melting and solidifi- 
cations. The order of magnitude of the shear modulus is a 
further source of complexity. The standard (zero-field) ex- 
pression for a body-centred cubic lattice (Puchs 1936) can 
be adapted, with equation (1), to give, 

1 1 w inl6 /vZ ^7-0.93 d1-6 „-3 rn\ 

/i = 1.1 X 10 Z^Z B12 erg cm . (9) 

However, the polar-cap gravitational constant is 3 2 x 10^* 
cm s"^ so that any density inversion may well induce a form 
of Ray leigh- Taylor instability. 

The final condensed-matter parameter that is impor- 
tant is the thermal conductivity parallel with B, which is 
extremely high (see Potekhin 1999). Thermal energy is de- 
posited at shower maxima a distance Zp below the surface 
z = 0, which is defined as that separating condensed matter 
from the atmosphere, and is then dissipated as blackbody 
radiation from the polar cap. The value of Zp in the high- 
density condensed matter of the neutron-star surface de- 
pends on shower energy owing to the Landau-Pomeranchuk- 
Migdal (LPM) effect but for the order of magnitude of the 
characteristic time we can assume —Zp ~ lOlr « 10~^ cm 
using the radiation lengths given by equation (2) or in Ta- 
ble 1 of Paper I. The characteristic time for shower energy 
input to produce a surface-temperature fluctuation is then, 

Cz^ 

teor^d = TTT^ « 10"® S, (10) 

2A|| 

in which typical values of the parameters are the specific 
heat C = 1.0 X 10^^ erg cm~^ and the longitudinal co- 
efficent of thermal conductivity Ay = 6 x 10^* erg cm"'^ 
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K^^. But for a surface temperature of 10^ K, the internal 
temperature gradient needed to conduct the radiated energy 
flux is extremely small, of the order of 10^ K cm~^, within 
the condensed matter at 2 < 0. In effect, heat is more easily 
conducted to greater depths than radiated from the surface. 
Consequently, the Green function G{z,Zp,t) giving the in- 
ternal temperature distribution must be almost independent 
of Zp and very close to that for zero temperature gradient at 
the surface z = 0. Thus a shower heat input H6{t) produces 
a temperature distribution within the condensed matter at 
« < 0, 

T(M) = //G(M)«(^^4^exp(^) (11) 

It is asymptotically oc t~^^^ so that the polar-cap tempera- 
ture arising from a sequence of heating events, each produc- 
ing a maximum temperature Tmax, has fluctuations away 
from T whose magnitude is a function of the time-scale con- 
cerned. Further discussion of this is deferred until our con- 
sideration of observed polar-cap blackbody temperature in 
Section 5. 



3 SHORT TIME-SCALE INSTABILITY 

Under the assumption that the case (ii) boundary condition 
is satisfied, the neutron-star surface has an atmosphere in 
local thermodynamic equilibrium (LTE) with approximate 
scale height Ia = {Z+l)kBT/Mg ~ 10~^ cm at temperature 
T = 10^ K, where Z and M are the mean ion charge and 
mass of the partially ionized atom. The expression for the 
chemical potential of an ideal Boltzmann gas gives the LTE 
atmospheric number density at z = 0, 

for atmospheric temperatures such that Na{0) <^ N. Its 
order of magnitude is Na{0) « 10^^ exp(-£c/A;BT) cm"^ 
for M = 56mp. We estimate the critical temperature Tc 
which must not be exceeded if the case (iii) boundary con- 
dition is to be valid by equating the flux of ions in such 
an LTE atmosphere that are incident on the neutron-star 
surface with the ion flux needed to give an open magneto- 
sphere current density equal to poc. This gives A''a(0) f» 10^'' 
cm~^ and ksTc = 0.025i?c. But the extent to which the 
atmosphere can be described as thin is very temperature- 
dependent. Thus for T = 2Tc and Z = 26, Na{Q) is of the 
order of 10^^ cm~^ and the whole atmosphere contains ions 
equivalent to about 10'^ Ir- At a density Na{0) > 10^'^ -10^" 
cm~^, the Boltzmann gas estimate of I a assumed here is un- 
reliable and it becomes necessary to allow for the interaction 
of protons with the ion atmosphere. 

If instabilities in plasma acceleration with time-scales as 
short as~ 10^* s are considered, the temperature at z = is 
constant apart from very small fluctuations, as shown by the 
Green function given at the end of the previous Section. Con- 
sequently, the temperature distribution and number density 
of ions in the atmosphere arc also constant in time and the 
atmosphere is in local thermodynamic equilibrium in the in- 
terval < z < zi, where zi is the top of the ion atmosphere, 
defined as the surface of last scattering. For time-scales of 



the order of Is, the temperature and depth of the atmo- 
sphere can change appreciably but there is no doubt that 
it is always in local thermodynamic equilibrium. The pro- 
cesses of proton formation by GDR decay followed by diffu- 
sion to the surface were described in Paper I, which assumed 
a very thin atmosphere. Therefore, we need to consider here, 
in more detail, diffusion in the atmosphere < z < zi. The 
number of protons is very many orders of magnitude smaller 
than that of ions so that the properties of the atmosphere, 
its scale height and equilibrium electric field given by the 
electrical neutrality condition, are determined solely by the 
latter. With the assumption of a single ion component we 
can determine the equilibrium electric field within the LTE 
atmosphere in the presence of a gravitational acceleration g 
and find that the proton potential energy is, 

^1 — mpQZ = ampQZ (12) 

for ions of charge Z and mass number A, and that this 
transports protons to 2 > «i from which region they are 
accelerated to relativistic energies in preference to the ions. 
Proton diffusion at a rate greater than is needed for a cur- 
rent density J'' = poc therefore cuts off ion acceleration and 
produces a thin electrically neutral proton atmosphere at 
z > zi. There appears to be no reason why this atmospheric 
structure should be disturbed by turbulent mixing. 

The proton average potential energy at 2 < is close 
to mpQZ but the jump bias, rripgas/kBT for ion separation 
as, is too small to be significant because nipgZp/kBT <§; 1. 
The atmospheric proton density at 2 = is necessarily some 
orders of magnitude smaller than the density at 2 < owing 
to the density discontinuity at the surface. Thus diffusion to 
the surface at 2 = is little changed by the presence of 
an atmosphere that is not necessarily very thin. Movement 
of the protons through the ion atmosphere is effected by 
the chemical potential gradient which is a function of ion 
number density. At high densities, as in the solid at 2 < 0, 
this is determined principally by entropy but as the density 
reduces, the potential energy given by equation (12) becomes 
the more important factor. In effect, the motion changes 
from a random walk to a drift velocity. The dividing ion 
density is given by the condition \r = Nj^ and is w 
4 X 10^^ cm~^, below which it is appropriate to define a 
drift velocity dotorminod by the local value of the moan free 
path \r derived from the total cross-section for Rutherford 
back-scattering, 

in which. 

The drift velocity is v ~ 2.6 cm s~^ at this density and 
T — lO^K. The consequence for a T = 2Tc atmosphere is 
that the diffusion time to 2 = 21 is increased, by some orders 
of magnitude, compared with the values assumed in Paper 
I which were for diffusion through the condensed matter 
only. Also, the change from random walk to drift velocity in 
the atmosphere has a significant effect on the distribution 
of diffusion times because it removes the long tail in the 
distribution that is a feature of random walks. 



6 P. B. Jones 



Some efTort has been made here to confirm that there 
is outward proton difi'usion because the process is of crucial 
importance. It was shown in Paper I that plasma accelera- 
tion is unstable, consisting of alternate bursts of proton and 
ion acceleration, with time-scale determined by the time for 
proton diffusion from z = Zp to z = zi. The polar-cap cur- 
rent density J is essentially time-independent for relativistic 
flow (as is the charge density) so that there is no fluctua- 
tion in polar-cap electric field other than that arising from 
the reverse flow of photo-electrons from accelerated ions (see 
Section 5). The instability is in the nature of the plasma ac- 
celerated and remains adequately described by the analysis 
given in Paper I which will not be repeated here. 

The maximum acceleration potential difference, given 
a current density J = poc, is almost entirely dependent on 
the Muslimov & Tsygan (1992) general-relativistic correc- 
tion. There is a separate contribution arising from ion iner- 
tia (Michel 1974) but it is important only at very low alti- 
tudes. Neglecting this, the max;imum potential difference is 
approximately 



271^ kB a 



7 X 10 



p2 



volts 



(15) 



and, unless restricted by pair creation at lower altitudes, is 
reached at an altitude z of the order of the neutron-star 
radius R which is roughly two orders larger than the polar- 
cap radius. We refer to Harding & Muslimov (2001, 2002) 
for a complete account. We assume here that for this current 
density, spontaneous pair creation by curvature radiation is 
not possible. Then the primary source of any positron com- 
ponent in J can only be the reverse flow of photo-electrons 
from accelerated partially-ionized atoms as discussed in Pa- 
per 1. But there is an essential difference here in that pro- 
tons in the very low density region at z ~ «i are almost 
completely ionized (here, we refer to Fig.l of Potekhin et al 
2006) so that both track length and energy flux of reverse- 
flow electrons from photo-ionized hydrogen atoms are negli- 
gibly small. Thus positron production in any interval of pro- 
ton acceleration is negligible. This is also true for ions of low 
atomic number Z ~ 4 — 5 which are completely ionized, but 
for higher Z, the reverse flux of inward accelerated electrons 
produces polar-cap ICS photons, as would outward accel- 
erated electrons in case (i). Pairs are produced by photons 
that are scattered to transverse momenta above the mag- 
netic conversion thresholds. The only difference is that the 
photons are inward moving. But even if spontaneous CR pair 
formation is not possible, positrons accelerated outward will 
produce CR pairs at higher altitudes, as in case (i), though 
superimposed on a flux of ions. 

In case (i), the reverse flux of positrons must form pro- 
tons which are not accelerated but form an atmosphere at 
z > whose equilibrium is defined by various diffusion pro- 
cesses perpendicular to the magnetic flux. Backward moving 
photons from the electron or positron showers arc a further 
source of pair creation in each of cases (i) - (iii) . These arise 
from the decay of residual nuclei following proton or neu- 
tron emission in GDR formation, and from (n, 7) reactions. 
Those photons that are not absorbed in the more dense part 
of the atmosphere can produce pairs if their transverse mo- 
menta exceed the threshold. The LPM effect, which in the 
high-density condensed matter at z < is significant for 
electrons of more than 10 GeV or photons of more than 2 



GeV, is less important in the atmosphere within which there 
will be some GDR formation. But we have been unable to 
make a satisfactory quantitative estimate of pair formation 
through this process. 



4 MEDIUM TIME-SCALE INSTABILITY 

Instability on time-scales some orders of magnitude longer 
than those of Section 3 is also of interest and can appear as 
a fiuctuation in the charge of nuclei reaching the polar-cap 
surface which we shall define here to be always at z = 
although it may move with respect to coordinates fixed at 
the centre of the star. Ions of initial charge Za move upward 
through the region of the shower maxima at Zp and, with 
nuclear charge reduced to Zb by GDR formation and decay, 
enter the atmosphere at « > 0. We wish to find if there 
are conditions under which the local average nuclear charge 
Zo{z) fails to be a time- independent function of position 
with limits Za ^ Zo{z) ^ Zb. The work of this Section 
replaces that of Section 4.1 in Paper I which contains an 
error. 

The longer time-scales here enable us to assume the 
proton and ion current densities and are the time av- 
erages of those described in the previous Section. There is, 
of course, a distribution of discrete nuclear charges in the 
shower region, but we shall work in terms of the average 
Z(z,t), the corresponding rmmber density N{z,t) given by 
equation (1), and the velocity v{z,t) with which nuclei ap- 
proach the polar-cap surface at 2 = 0. Proton formation by 
GDR decay occurs predominantly in the very late stages of 
shower development, as explained in Section 1, and so we 
shall assume that it is confined within limits Za and Zb and 
is defined by the normalized function gp{z), 



gp{z)dz ■■ 



(16) 



The physical basis for our study of the system is that the 
total rmmbors of nuclei and of protons (bound or unbound) 
are conserved. Because both neutrons and protons are pro- 
duced in GDR decay, we can make the approximation of 
neglecting the effect of /3-transitions. Therefore, the conti- 
nuity equations are. 



dN d{Nv) _ 
dt ^ dz ~ ' 



(17) 



and with neglect of the relatively short proton diffusion time 
so that within medium time-scales all protons produced in 
the shower are assumed to be promptly accelerated from the 

atmosphere at 2; « 21, 



^(^^)+^=-,,(2)J^(t). 



dt 



dz 



(18) 



Because the reverse electron flux from photo-electric ioniza- 
tion is the source of the showers, we shall find it convenient 
here, as in Paper 1, to introduce the parameter K which 
is a function of the atmospheric rmclear charge and is the 
number of protons produced per unit positive nuclear charge 
accelerated. With the representation of equation (1) in the 
form N = CZ'' , equations (17) and (18) can be combined 
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to give, 

C r dz(z-<^+vZ-'^)=-J-{t). (19) 
It has the obvious time-independent solution, 



Za — Zb = KZf,, 



(20) 



and the time-independent velocity of nuclei is such that 
Z^{z)vo{z) is independent of z. 

A natural fluctuation away from Zq would be of the 

form Z{z, t) = Zq + 5Z{z, t) with, 

5Z = 'n{Za-ZQ{z))e"'\ (21) 

giving similar fluctuations in J'', N and v, in which r) is 
infinitesimal and independent of z and t. Prom equation (17) 
we then have, 

5{Z''v) = — iwr;x(a)e"^*. 



where. 



X(«)=7 / dz' {Za - Zo) Z^-\ 

^ Zn. 



(22) 



(23) 



Substitution into equation (19) with the retention of terms 

of first order in r\ gives, 

iwTje'"* / ' dz (zl {Za - Zo) - xiz)- 



^dZo\ 
dz J 



Z^vo (Zt - Za) r,e'^* = -^SJ^{t) 



(24) 



The current densities, averaged over short time-scales, arc 
J"^ = NZv{0,t) and, from the definition of K, J" = KJ\ 
The fluctuation away from the time-independent solution is 
5J — SJ^ + SJ^ = 0, which wc assume to be maintained by 
the boundary conditions on <1>. To express 5J^ in terms of 
5Z we represent the Z-dependence of K in the vicinity of Zi 
as a power law K = KoZ'^(Q, t). We find, after eliminating 
the velocity fluctuation 5v{Q,t), that 



yr = yZ^^v,^^5Z(G,t). 



(25) 



The relationship with equation (24) is established by noting 
that, 



<5Z(0, t) = 5Z{zb, t-T) = r}{Za- Zk) e 



iw(t-T) 



(26) 



where r is the time interval of rmclcar movement from z — Zf, 
to « = 0. From equations (24) - (26), the equation whose root 
u) we require can be expressed as, 



Z^vo{zb) 



f 



dzZ^-' ((1 + 7)Zo - -yZt) {Za - Zo) 



= -{Za - Zb) - 



vK 



{Za - Zb)e- 



(27) 



K + V 

in which the function x{z) has been removed by partial in- 
tegration. 

The only assumptions we need malce about the depth 
distribution of proton formation in the late stage of shower 
development are that it is small outside the interval Za < z < 

Zb and that z^ — Za is smaller than \zb \ though not necessarily 
so by as much as an order of magnitude. A suitable function 
would be. 



I W 2 . 2 ( 7r(0 — Za) 

gp{z)) = sm ' 



Zb — Za 



Zb — Za 



(28) 



in terms of which we could express Zo as an explicit function 
of z, which would be needed to obtain numerical values for 
the root u! = uji + iw2 of equation (27). However, we are 
primarily concerned here not with growth rates but with 
the boundary between stability and instability and so shall 
not proceed with numerical solutions for w. Fortunately, it is 
possible to obtain a sufficient condition for the existence of 
instability independently of the form of Qp. Provided —1 < 
7 < 0, which is clearly satisfled, we can see directly from 
the values of the integrand in equation (27) at the limits 
and from its lack of an extremum that. 



Zn. 



dz 



< {Zb - Za){Za - Zb), 



(29) 



and therefore that equation (27) can be replaced by the in- 
equality. 



Zb 



> -1 - 



vK 



\Zb\ I K+1 
Thusa;i,2 satisfy the inequalities. 



(30) 



ui2r 



Zb - Za 
\Zb\ 

Zb - Zg 



> 



K + 1 

I 

K+1 



e ^ smwiT 



< 1 + — — -e coswiT. 



(31) 



From the first of these, wc can see that the real part of 
ui must satisfy cjir > 7r/2, provided the proton formation 
region in the shower is sufficiently compact that. 



vK 

kT-i' 



> 



Zb ■ 



Zb 



(32) 



With greater compactness, loit — ^ n, as would be expected. 
The second inequality gives the condition for 012 < 0, that 
is, fluctuation growth. At the threshold where \uj2t\ <C 1, 
and for Zb — Za \zb\, it is simply vK > K + 1. 

The Z-dependence of K was represented in Paper I as 
an approximate power law, K — KqZ" with v = 0.85 for 
Z ^ 10 but with the reservation that this would certainly 
become invalid for an atmosphere of ions with Z ~ 5 which 
would be almost completely ionized. Therefore, we should 
anticipate quite large values of v at small Z, sufficient to 
give the unstable behaviour found here. Stability clearly de- 
pends on the energy flux from photo-ionizatiou. Low fluxes 
and moderate rates of proton formation in the GDR region 
of the showers, such that the nuclear charge inferred from 
equation (20) is Zb ~ 10, allow a stable time-independent 
progression of nuclear charge as a function of depth below 
the polar-cap surface as given by Zq{z). But larger values 
of K and smaller Zb lead to instability. (The value of Za is 
unimportant, provided it is not small, and for order of mag- 
nitude estimations in this paper we have assumed Za = 26 
unless otherwise stated.) 

Unfortunately, our analysis of the instability of nuclear 
movement to the polar-cap surface is limited to small fluc- 
tuations and does not extend to large amplitudes. But it is 
not difficult to see the form it would take. An atmosphere of 
high-Z ions at time t produces a high reverse-electron energy 
flux which creates a layer of low-Z^ ions in the shower max- 
imum region Za < z < Zb- These flow toward the surface at 
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= and form a \ow-Z atmosphere at time t + r which pro- 
duces a low reverse-electron energy flux and correspondingly 
large values of Zi,. The ions accelerated alternate between 
high and low-Z values. The basic unit of time is given by 
the time tri for the emission at the Goldreich-Julian current 
density of one radiation length of ions, 

tri = 2.1 X W^Z'^B^^i-Pscc'iJj) 

(ln(l2Z^/'^B-2^''^^ ' s. (33) 

The high-Z intervals are subject to short time-scale instabil- 
ity as described in Section 3, but low-Z intervals have little 

or no reverse-electron flux and therefore no possibility of 
significant positron acceleration and electron-positron pair 
production. 

The instability outlined here is of quite simple form, but 
there are complicating factors that have been mentioned ear- 
lier in the Section 2 consideration of polar-cap parameters. 
Evaluations of the melting temperature discussed following 
equation (8) show that there is every possibility that the 
condensed matter state below the atmosphere may be liq- 
uid, or a solid close to melting with a high self-diffusion rate. 
This would have no effect on short time-scale instability but 
could complicate the behaviour of the system over medium 
time-scales. The melting temperature is Z-dependent so that 
a density inversion is possible, the upper layer of highor- 
Z being either liquid or solid. In the liquid case, we must 
anticipate Rayleigh- Taylor instability, which may also exist 
in the solid case because its shear modulus (see equation 
9) may not be adequate to maintain mechanical stability. 
All these processes are occurring at depths \zp\ ~ 10~^ cm 
but over a polar cap whose radius is of the order of 10'' 
cm. Consequently, a further complication is that diflferent 
polar-cap areas are unlikely to be in phase with each other. 
These problems are considered further, though necessarily 
in a qualitative way, in Section 5. 



5 NULLS SUBPULSE DRIFT AND 
POLAR-CAP COHERENCE 

The time-dependent phenomena considered in Sections 3 
and 4 are local and one-dimensional because both shower 
depth Zp and atmospheric scale height are very many or- 
ders of magnitude smaller than the polar-cap radius. This 
introduces the question of whether or not there is coherence 
over the whole polar-cap area. Both instabilities are depen- 
dent on the parameter K which is a function of the surface 
nuclear charge Z{0, t) and to a lesser extent of surface tem- 
perature and acceleration field. For this reason alone, we see 
no case for assuming complete coherence and anticipate that 
the polar cap we describe has zones of proton and ion emis- 
sion which cannot be stationary, the total areas of each being 
determined, approximately, by the average value of K. For 
neutron stars that are unable to support spontaneous pair 
production by curvature radiation, the proton zones have 
no reverse electron flux, do not support pair production, and 
hence merely produce an accelerated one-component plasma 
as described in Section 3. But the ion zones will support ICS 
pair production and so appear as moving sources within the 
polar cap. In this Section, we shall attempt to compare pos- 



sible forms this motion might take with some of the observed 
phenomena in radio pulsars. 

The distinction between the average pulse profile and 
individual sub-pulses within it was noted almost immedi- 
ately following the discovery of pulsars (sec, for example, 
Smith 1977). The amplitude and form of successive sub- 
pulses can vary in times of the order of the rotation pe- 
riod and in some pulsars, observed with higher resolution, 
sub-pulses have micro-structure of 10^^ — 10^'^ s time-scale. 
There are also more organized phenomena, and for recent ex- 
tensive surveys of these we refer to Wang et al (2007) in the 
case of pulse nulls and to Weltevrede, Edwards & Stappers 
(2006) for sub-pulse drift. There is a move toward a consen- 
sus (Lyne ct al 2010) that quasi-periodic switching between 
magnetosphoric states with different spin-down rates is the 
basis of mode-changing. These authors even suggest that it 
is the source of a large component of pulsar timing noise. 
However, it is also true that the subpulsc characteristics ob- 
served in a small number of pulsars are quite singular, but 
here discussion is confined to the general features of sub- 
pulses. 

Given the properties of the medium time-scale instabil- 
ity described in Section 4, it is natural to associate with nulls 
those intervals in which high-Z ion zones either do not exist 
or are confined to areas of the polar cap firom which radia^- 
tion produced by the plasma is not visible to the observer. 
For sufficiently low values of Z^, ions are accelerated from 
the surface completely ionized so that there is no reverse- 
electron flux and hence no pair creation and radio emission. 
The current density is little changed in the open sector of 
the magnctosphoro but the absence of pair creation means 
that the particle content near the light cylinder is quite dif- 
ferent as are the components of the momentum density or 
stress tensor on any spherical surface in this region centred 
on the star. It is therefore not surprising that the spin-down 
torque is reduced during the interval of a null, as has been 
observed in PSR B1931-I-24 (Kramer et al 2006). Neutron 
stars with small K such that k. 10 are likely to have 
V < 1 and so will have a steady-state progressively reduc- 
ing nuclear charge Zo{z), no medium time-scale instablity, 
and no long-term nulls. But the questions about melting and 
mechanical stability of the polar cap which were mentioned 
at the end of Section 4 remain and the whole system is so 
complex and difficult to describe in physical terms that we 
are unable to give useful quantitative predictions. However, 
incomplete nulls, having a low but detectable level of emis- 
sion should be observed. It is also unsurprising that nulls are 
observed to be not completely random (Redman & Rankin 
2009). 

To some extent and unfortunately, the same remarks 
have to be made about the effects of short time-scale insta- 
bility. For temperatures T > 2Tc, the proton diffusion time 
is much longer than in the case of the very thin atmosphere 
assumed in Paper I but remains difficult to calculate with 
complete confidence. In order to describe the polar cap it 
is necessary to introduce coordinates u{z) on a surface per- 
pendicular to B. As in Paper I, the proton current density 
at any point u(0) can be related to the ion current density 
through the definition of K, 

J''iu,t) + J''iu,t) = K [ dt'fpit-t')J^u,t'), (34) 
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in which fp is the normaUzed proton difTusion-time distribu- 
tion. Without the J'' term, whose significance is described 
below, this would be a homogeneous Volterra equation of the 
second kind having no non-zero square-integrable solution. 
(The timc-dcpcndcnco of K arising from the temperature de- 
pendence of the LTE ionic charge Z is neglected here.) The 
approximate expression for fp given in Paper I (equation 
23) assumed a random-walk diffusion through the condensed 
matter at z < and so is not valid for an atmosphere with 
properties given by equations (13) and (14). For T > 2Tc, 
most of the diffusion time is in the drift-velocity phase, in 
which case the time distribution would be more suitably ap- 
proximated by a normalized gaussian function centred at 
t — Tp or, in the limit, by fp{t — t') = 5{t — t' — Tp), where 
Tp is derived from equation (13) of Section 3. The quantity 
JP = within intervals for which < 1 and at other times, 
when J'' = J, represents the storage of excess protons in the 
atmosphere at 2; > The total current density J is fixed by 
the boundary conditions and at 2 = differs little from the 
Goldreich- Julian value. Thus = J and = until the 
instant at which the atmosphere is exhausted and falls to 
some residual value < J. Then ion flow recommences al- 
most immediately (in a time much shorter than the growth 
time for spontaneous curvature radiation pairs) and contin- 
ues until proton diffusion grows sufficiently to re-form the 
proton atmosphere. The durations of the time intervals for 
ion and proton emission are those for electron-positron pair 
creation or otherwise, and labelled Tee and Tgap are given by. 




-t')r{t') 



Tgap = -J J dt'r{t') (35) 

Therefore Tee and Tgap both depend on the smallness of fp 
for small values oi t — t' . In the (5-function limit for fp, wc 
have Tgap = KTee- Estimates of the parameter K = KoZ" 
were given in Paper I, but are repeated here, 

Ko = 2.8(Z.^)-°-^«SiTTe-^ °, B12 > 1, 

Ko = 1.6{Z,mr°''^Ti^°, Bi2<l, (36) 

with f = 0.85 for 10 < Z < 26. In these expressions, T is not 
the local surface temperature but is an average for radiation 
emitted over the whole polar-cap, and B is the actual polar- 
cap magnetic flux density. For B12 S> 1, the values of Ko 
given need to be modifled, though not greatly, to allow for 
the high-B values of Wp contained in Table 1. 

Guided by the estimates contained in equations (6) and 
(7), we anticipate that except for a very small number of neu- 
tron stars, the surface temperature is at all times T > Te so 
that the case (ii) boundary condition is always maintained 
and ion emission is never temperature-limited. It is limited 
instead by the fact that the proton atmosphere forms above 
the ions and the protons are preferentially accelerated, as 
described in Section 3. The only effect of the temperature 
variations described by equation (11) that occur as a result 
of alternating proton and ion emission is to change the den- 
sity of the LTE ion atmosphere. It is also worth observing 
that the local nature of equation (34) is not disturbed by 
the presence of E x B drift above the polar cap. Although 
this slightly displaces a reverse electron shower relative to 



the point from which the ion was accelerated, it has no effect 
on ion emission, which is not temperature-dependent. 

In view of the time-variation described by equations 
(34) and (35), is it possible to imagine organized rather than 
chaotic motion of ion zones on the polar cap? An example 
of chaotic motion would be the existence of very many small 
zones without obvious organization. Two simple organized 
cases would be motion along a slot of constant rotational lat- 
itude and circular motion at constant m(0). With the Desh- 
pande & Rankin (1999) analysis of PSR B0943+10 in mind, 
we consider circular motion. Equations (34) and (35) are 
local in u and contain no information that can determine 
an ion zone movement velocity u. The quantities that are 
essentially constant are K and, apart from the effect of vary- 
ing LTE atmospheric temperature, Tp. Thus the circulation 
time for n ion zones, is P3 = nPz, where P3 is the band 
separation in the usual notation by which subpulse drift is 
described, and here is given by 

P3 = Tgap + Tee = (K + l)Tp (37) 

in the 5- function limit of the diffusion function fp. In this 
system, the velocity ii is determined by n for fixed values 

of K and Tp. For the same reason, drift direction is also un- 
specified and there is nothing to preclude a change in n or a 
reversal following some disturbance to the polar-cap surface 
such as might be a consequence of the kind of mechanical 
instability briefly described in Section 4. The drift time is 
determined principally by diffusion through the more dense 
layers of the atmosphere. From equations (13) and (14) it is. 




In Section 3 we observed that atmospheric density at 2 = is 
an exponential function of T > Tc so that its value is essen- 
tially unpredictable. Therefore, it is likely that the density 
discontinuity at z = is small and that the reverse-electron 
showers may be contained entirely within the atmosphere. In 
this case, the lower integration limit in equation (38) must 
be replaced by Zsm > 0. The diffusion time is then not di- 
rectly a function of B, but is dependent on surface gravity 
and on Z. But its distribution for all pulsars should be com- 
pact. The distribution in the values of K, given by equation 
(36), is probably the more important source of the observed 
spread in the values of P3 . Detailed calculation of Tp has not 
been attempted here. In particular, our estimates of Xr and 
Ia are subject to some uncertainty. 

The fact that P3 is constant for a particular pulsar 
whereas the circulation time P3 is dependent on n could in 
principle allow comparison with E x B drift-velocity polar- 
cap models in which P3 is constant. It is also worth consid- 
ering the effect of mechanical instability in allowing leakage 
of protons (or law-Z ions) to limited areas of the surface as 
described at the end of Section 4. The excess protons form 
a localized atmosphere of greater depth than in surrounding 
areas, so that ion zones (with consequent pair creation) do 
not form there until it is exhausted. The observer of radio 
emission probably sees plasma from no more than a strip of 
polar-cap area at roughly constant rotational latitude so it is 
possible that intervals of null emission would be be observed 
at fixed longitude in a sequence of subpulse bands. 

In a survey of 187 pulsars, selected only by signal-to- 
noise ratio, Weltevrede et al (2006) found that reversals do 
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occur and that roughly equal numbers of pulsars have sub- 
pulse drift to smaller or greater longitudes. They also com- 
ment that subpulse drift is so common a phenomenon that 
it cannot depend on pulsar parameters having extraordi- 
nary values. The band separation P;^ is independent of radio 
frequency, and is not correlated with period P, age P/2P, 
or with the inferred dipole field Bd- Measured values of P3 
given for 77 pulsars in Table 2 of Weltevrede et al are mostly 
contained within a single order of magnitude, 1 < P3 < 10 
s although the distribution has a small tail extending to 
~ 20 s. This is quite compact (for a neutron-star parame- 
ter) and is not inconsistent with the interpretation of sub- 
pulse drift given here and with equation (37). The range of 
P3, with estimates of K found from equations (36), indicate 
Tp ~ 10"^ - 10° s and Na{Q) = 10^^ - lO^'* cm'^', values 
that are by no means unreasonable. 

It is necessary to compare our polar-cap model with 
two other sets of pulsar observations. Given that protons 
are the major fraction of particles accelerated and produce 
no electron-positron pairs, it is worthwhile considering the 
radio luminosity L^, . The order of magnitude of this can be 
expressed as, 



L^AOAf « — — e i 

ceP^ 



1.4 X l()"'"lj,fi2 
p2 



(39) 



in terms of the neutron-stax radius R and the inferred dipole 
field Bd- Prom this we can estimate the energy e radi- 
ated into solid angle Ail within bandwidth /S.i> per unit 
charge (baryonic or leptonic) accelerated at the polar cap. 
Using the 400 MHz luminosities hsted in the ATNF catar 
logue (Manchester et al 2005) and a bandwidth Ai/ = 400 
MHz, we find by evaluating equation (37) for a small sam- 
ple (B0826-34, B0834+06, B0943+10, B0950+08, B1055-52, 
B1133-I-16, B1929-t-10) that typical values arc in the inter- 
val € = (40 — 400)(5n MeV. Even though the emission solid 
angle is likely to be as small as Af2 ~ 10^"^ — 10^^ stcrad, 
there rrmst be some concern here because the number of ICS 
pairs produced per primary accelerated positron (Harding & 
Muslimov 2002) is not large so that either the conversion of 
electron-positron energy to coherent radio emission is effi- 
cient or other plasma components are involved, as in the 
paper by Cheng & Rudcrman (1980). 

A second set of observations that arc relevant are those 
of polar-cap blackbody X-ray luminosities. The blackbody 
X-ray emission expected from a polar-cap of ion and proton 
zones can be found from equation (11). The reverse-electron 
flux from photo-ionization heats an ion zone at a rate Hq 
within a time interval < t < Tp. An estimate of this can 
be obtained directly from the mean electron energy per unit 
nuclear charge accelerated, given by equation (29) of paper 
I, and is, 

Ho = 6.0 X 10^*Z° **^(0)PiVT6"'P"\ (40) 

in units of erg cm~^ s~^. In this expression, T is not the local 
surface temperature but is an average for radiation emitted 

over the whole polar cap. The surface temperature derived 
from equation (11) with neglect of radiative loss is 

2Ho 



T{t) = 



(7rCA||)i/2 



1/2 



< t < Tj,, 



(41) 



and increases very rapidly, so that the limiting temperature, 
T06 w 0.6Z°■l^(0)B?2^p-°■^ (42) 



derived from Hq is reached ai t <^Tp. The cooling ai t > Tp 
is also rapid so that the major part of the X-ray luminosity 
is that of a black body at To and of area approximately 
equal to the canonical dipole-field area 2'k'^P? jcP divided 
by 7f + 1. 

There have been many attempts to measure the polar- 
cap blackbody temperature and source area of a sub- 
set of radio pulsars. We refer to Zavlin & Pavlov (2004) 
for B0950+08 ; De Luca et al (2005) for B0656-hl4 and 
B1055-52; Tepcdelcrflioglu & Ogclman (2005) for B0628-28; 
Zhang, Sanwal & Pavlov (2005) for B0943-M0; Kargaltsev, 
Pavlov & Garmire (2006) for B1133-M6; Gil et al (2008) 
for B0834+06; Misanovic, Pavlov & Garmire (2008) for 
B1929+10. Five pulsars (B0628-28, B0834-f06, B0943+10, 
B1133+16, B1929+10) have source areas one or two orders 
of magnitude smaller than the canonical area, but unfortu- 
nately, a systematic comparison with equations (40)- (42) is 
not possible because in most instances the authors are able 
to say only that the observed X-ray spectrum is consistent 
with the stated temperature and source area. The quoted 
source temperatures are ~ 3 x 10® K and are larger than 
those predicted by equation (42) unless it is assumed that 
B ^ Bd- They arc also uncomfortably large in the context 
of E X B drift-velocity polar-cap models such as that devel- 
oped by Gil, Mclikidze & Geppert (2003) although it must 
be conceded that this model allows an interesting test of its 
validity (see Gil, Melikidze & Zhang 2007). However, the 
case (iii) surface electric-field boundary condition on which 
these models rely could be maintained only for ion cohesive 
energies of ~ 10 keV which in turn would imply actual fields 
two orders of magnitude or more larger than the ATNF cat- 
alogue dipole field. 



6 CONCLUSIONS 

This paper is a continuation of a previous study (Paper I) of 
isolated neutron stars with positive corotational charge den- 
sity and surface electric-field boundary condition E • B = 
at the polar cap. The reverse flux of electrons arising 
from photo-ionization of accelerated ions is incident on the 
neutron-star surface and produces protons through forma- 
tion and decay of the giant dipole resonance in the later 
stages of electromagnetic shower development. Protons are 
the major component of the accelerated plasma but we flnd 
that a time-independent composition of ions, protons and 
positrons is usually unstable. The consequences of these phe- 
nomena should be observable in pulsars unless, of course, 
rendered nugatory by some factor not properly taken into 
account here. But we emphasize that there are almost cer- 
tainly other sources of instability contributing to the com- 
plex behaviour observed in radio pulsars that may be present 
in each of cases (i) - (iii). 

There are two instabilities which result in transitions 
between states of different plasma composition above local- 
ized areas on the polar cap and it is suggested here that these 
are the basis for the commonly observed phenomena of nulls 
and of subpulse drift. The beisic units of time are Tp ^ 
s for the short time-scale instability and for the medium 
time-scale, tri ^ 10^ — 10^ s. These instabilities are not 
primarily electromagnetic in origin and will not be present 
in neutron stars with negative corotational charge density 
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and electron acceleration. Micropulscs of ~ 10^^ — 10^'^ s 
duration are also present in some pulsar emissions but arc 
unlikely to be associated with the instabilities we consider 
here. There appears to be no reason why both signs of coro- 
tational charge density should not be present in the isolated 
neutron-star population and it would be interesting to see if 
there is any observational evidence for an associated division 
of the ATNF catalogue pulsars. 

The transitions are principally in plasma composition 
and the current density at the neutron star surface remains 
close to J = poc. The electric potential ^{z,u{z)) in the 
open magnetosphere above the polar cap is also almost un- 
changed. These transitions are, in principle, observable be- 
cause there are no processes by which electron-positron pair 
creation can occur in a plasma containing only protons, 
whereas it does occur in the ion plasma. Thus there will be 
no coherent radio emission in the proton phase if we assume, 
as is usual, that a pair plasma is a necessary constituent for 
it. 

Both the medium time-scale instability and the mechan- 
ical instability of the neutron-stax surface mentioned in Sec- 
tion 4 provide a basis for pulse nulls, but unfortunately, not 

one that is suitable for quantitative prediction. The short 
time-scale instability is more interesting. We show that at 
any instant, the polar cap is divided into moving ion and 
proton zones and consider organized motion of these, specif- 
ically the circular movement of compact ion zones around 
the magnetic pole. We make the hypothesis that such mo- 
tions can occur (following recent work by Deshpande & 
Rankin 1999) but have not attempted to show that they 
have long-term stability, to the extent indicated by observa- 
tion, against decay to a chaotic state. This is of some interest 
because it shows that features usually considered to require 
the surface-field boundary condition E • B ^ are possible 
under the E • B = condition. Actual polar magnetic flux 
densities exceeding 10^'* G arc necessary for the former con- 
dition and it must be doubtful that such fields, two orders 
of magnitude greater than the inferred (catalogue) dipole 
field, are likely to exist in the very large number of pulsars 
observed by Weltevrede et al (2006) to show subpulse drift. 
Our view of subpulse drift is also interesting in that the 
motion is not an E x B drift velocity as in the model of 
Ruderman & Sutherland. The fixed parameters for a given 
neutron star are the band separation Ps, and K which is the 
number of protons produced per unit ion charge accelerated. 
This latter parameter approximately determines the ratio of 
the total areas of proton and ion zones on the polar cap at 
any instant. Thus the circulation time P3 = nPs depends on 
n, the number of ion zones, and is not necessarily constant 
as it would be in the E x B drift velocity model. The band 
separation P3 is dependent on K and on the diffusion time 
Tp and so is independent of rotation period P and almost 
independent of B. 

We have assumed here, following the pair formation 
calculations of Hibschman & Arons (2001) and Harding & 
Muslimov (2002), that spontaneous pair creation by curva- 
ture radiation is not possible in most isolated neutron stars. 
Those with a sufficiently large value of |$maa;| oc Bd/P^ 
to support spontaneous CR pair creation, predominantly 
young high-field pulsars, are not expected to show the in- 
stabilities we have considered here. These authors actually 
considered case (i) in which the basis of plasma formation 



is an electron current density close to J = poc at the sur- 
face accelerated to energies sufficient for ICS or spontaneous 
CR pair formation. In case (ii), the protons and ions form 
the basic current component with |<I>maa;| broadly the same 
as in case (i) apart from the ion inertia term which is very 
small for electrons. Thus the condition for the spontaneous 
growth of CR pair formation is very similar in cases (i) and 
(ii). Although the reverse electron fltix arising from CR pair 
formation in case (ii) is probably quite a small fraction of 
poc, it is likely that, for the large values of |$maa;| that arc 
necessary, there will be a proton atmosphere which is never 
exhausted. In this instance, steady-state plasma formation 
and low-altitude acceleration appear probable, though sub- 
ject to possible instabilities at those higher altitudes where 
the coherent radiation is formed. The same statement can 
be made about all case (i) pulsars whether or not sponta- 
neous CR pair formation is supported. Microstructure with 
10~^ — 10~^ s timcscalos is likely to be a consequence of such 
higher-altitude instabilities. The reason for this assumption 
is that the case (i) boundary condition on $ needs only an 
electron density on the surface separating open from closed 
magnetospheres. The small electron cohesive energy means 
that there is no reason why this should not be maintained 
at all times, as is also true for the Goldreich- Julian elec- 
tron current density. For this boundary condition, there is 
no obvious way in which non-electromagnetic time-constants 
can influence plasma formation and acceleration at altitudes 
less than ~ lO" cm above the polar cap. It is possible, per- 
haps, to assign either this boundary condition or case (ii) 
with spontaneous CR pair formation to those pulsars in the 
systematic study of Weltevrede et al (2006), approximately 
one half of the total in number, that do not exhibit subpulse 
drift. 
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